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Abstract 

In a recent paper Jones and Mateo used operator techniques to show that the non-Hermitian 
PT-symmetric wrong-sign quartic Hamiltonian H = ^p 2 — gx 4 has the same spectrum as the 
conventional Hermitian Hamiltonian H = \p 2 + 4gx 4 — ^flgx. Here, this equivalence is demon- 
strated very simply by means of differential-equation techniques and, more importantly, by means 
of functional-integration techniques. It is shown that the linear term in the Hermitian Hamiltonian 
is anomalous; that is, this linear term has no classical analog. The anomaly arises because of the 
broken parity symmetry of the original non-Hermitian PT-symmetric Hamiltonian. This anomaly 
in the Hermitian form of a PT-symmetric quartic Hamiltonian is unchanged if a harmonic term 
is introduced into H. When there is a harmonic term, an immediate physical consequence of the 
anomaly is the appearance of bound states; if there were no anomaly term, there would be no 
bound states. Possible extensions of this work to — 4> quantum field theory in higher-dimensional 
space-time are discussed. 
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I. INTRODUCTION 



In this paper we consider the quantum system described by the Hamiltonian 

H = ±p<- gX \ (1) 

where g is real and positive. The Hamiltonian (pQ) is of particular interest because the 
corresponding — <f) 4 quantum field theory might be a good model for describing the dynamics 
of the Higgs sector of the standard model. This is because the — 4 theory is asymptotically 
free and thus nontrivial 0, 0, HI- Furthermore, the one-point Green's function (O|0|O) is 
nonvanishing ji]]. 

While the — x A potential in (Q) appears to be unbounded below and the Hamiltonian H 
appears to be symmetric under parity reflection, neither of these properties holds because, 
as we will explain shortly, the eigenfunctions of the Hamiltonian are required to vanish 
exponentially as |x| — > oo in a pair of Stokes' wedges in the complex-x plane. These wedges 
do not include the real-x axis. 

The Hamiltonian (JI} is not Hermitian in the conventional Dirac sense, where Hermitian 
conjugation is defined as combined matrix transposition and complex conjugation. Never- 
theless, the eigenvalues E n are all real, positive, and discrete. This is because H possesses 
an unbroken VT symmetry j^, 0, 0], which means that H and its eigenstates ip n (x) are in- 
variant under space-time reflection. Here, V denotes spatial reflection p — > —p and x — > —x, 
and T denotes time reversal p —>■ —p, x — > x, and i — ► —i. 

Even though H is not Dirac Hermitian, it is possible to construct a state space having a 
positive inner product with respect to which H is Hermitian and in this Hilbert space, time 
evolution is unitary (probability conserving) [f| . This positive inner product involves CVT- 
conjugation, where C is a linear operator whose square is unity and which commutes with 
the Hamiltonian. (These mathematical properties are reminiscent of the charge conjugation 
operator in particle physics.) Given the operator C, we can construct the positive operator 
e Q _ Q<p^ w hich can in turn be used to construct by means of a similarity transformation 
an equivalent Hamiltonian H = e~®/ 2 He®/ 2 , where H is Dirac Hermitian |8[. 

Many PT-symmetric quantum-mechanical Hamiltonians have been studied in the recent 
literature However, the Hamiltonian (Q) is especially interesting because, unlike the VT- 
symmetric Hamiltonian ^p 2 +ix 3 , for example, the boundary conditions on the eigenfunctions 
are not imposed on the real-x axis. Rather, these boundary conditions hold in the interiors 
of wedges in the complex-x plane. To identify these wedges we use a WKB approximation to 
determine the possible asymptotic behaviors of the eigenfunctions. For large \x\ the possible 
exponential behaviors of the solutions to the Schrodinger equation are given by 

^ n (x) ~ e^^" 3 / 3 (|x| -> oo). (2) 

This result shows that the eigenfunctions are purely oscillatory along six lines separated by 
60° angles. In particular, the wave functions are oscillatory on the positive and negative 
real- a; axes. These six lines are the boundaries of the six 60° Stokes' wedges in which the 
solutions in (J2J) grow or decay exponentially. 

The eigenfunctions for the Hamiltonian JI} are required to decay exponentially in the 
interiors of a pair of Stokes' wedges in the lower-half x-plane. These wedges, which are 
symmetrically placed with respect to the imaginary axis, lie below the positive and negative 
real- a; axes with the upper edges of the wedges lying on the real axis. 
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Note that under space reflection x — *> —x, the original two-wedge domain changes to the 
interior of a different pair of 60° wedges that lie in the upper-half plane. Therefore, while 
H in may appear to be parity symmetric, it obeys parity- violating boundary conditions. 
One can also understand the parity violation of H in (JTJ in a different way. In Ref. 
it is shown that the real spectrum of H in (|1J) can be obtained without having to impose 
boundary conditions in the complex plane. The procedure is simply to require that the — x 4 
potential be reflectionless. Thus, an incoming plane wave from x = — oo propagates past the 
potential and becomes an outgoing plane wave at x — +00. This configuration is not parity 
invariant because under parity reflection, the incoming plane wave at x = — 00 becomes an 
outgoing plane wave and the outgoing plane wave at x = 00 becomes an incoming plane 
wave. In short, the right-going flow of probability current becomes a left-going flow. (Of 
course, this configuration is invariant under combined V and T reflection.) 

This violation of parity symmetry occurs only at the quantum level. At the classical 
level, the equations of motion are clearly parity symmetric: 

x(t) = —p(t), p{t)=4gx 3 (t). (3) 

In general, it is difficult to construct the operator C, and in the past for nontrivial models 
this operator has only been determined perturbatively [11]. However, recently Jones and 
Mateo used perturbative operator methods to construct C in closed form for the Hamiltonian 
in (0) They then found the equivalent Hermitian Hamiltonian H, whose potential has 
a positive quartic term and also a linear term. (Using Rayleigh- Schrodinger perturbation 
theory, Buslaev and Grecchi had already discovered this equivalent Hamiltonian much earlier 
Q.) 

The purpose of this paper is to examine the connection between the PT-symmetric non- 
Hermitian Hamiltonian (JJJ and the equivalent Hermitian Hamiltonian H, whose potential 
has a positive quartic term and a linear term. In Sec. Owe demonstrate the equivalence be- 
tween these two Hamiltonians simply and directly by transforming the Schrodinger equation 
for the former Hamiltonian into the Schrodinger equation for the latter Hamiltonian. We 
show that the linear term in H is proportional to h and is a therefore quantum anomaly. This 
anomaly arises because the boundary conditions on the PT-symmetric Hamiltonian violate 
parity. As explained above, this parity violation is not a feature of the classical equations of 
motion. In Sec. lllll we generalize the Hamiltonian H in ([T} to include a harmonic (quadratic) 
term and show that the equivalent Hermitian Hamiltonian has an anomaly of exactly the 
same form. Furthermore, we show that an immediate consequence of the anomaly is the 
appearance of bound states. 

The Schrodinger-equation approach of Sec. |H] does not readily generalize to quantum 
field theory, so in Sec. |IV| we demonstrate the equivalence between H and H by using 
path-integration techniques in which we treat H as defining a quantum field theory in one- 
dimensional space-time. We conclude in Sec. |V] by indicating how path-integral techniques 
might be extended and used to identify the equivalent Hermitian Lagrangian for a —(f) 4 
quantum field theory in higher space-time dimensions. 



3 



FIG. 1: Stokes' wedges in the lower-half complex- a; plane for the Schrodinger equation (jlj) arising 
from the Hamiltonian H in ((J). The eigenfunctions of H decay exponentially as |x| — ► oo inside 
these wedges. Also shown is the contour in (JSJ. 



II. DIFFERENTIAL-EQUATION DERIVATION OF THE ANOMALY 



In this section we use straightforward differential-equation techniques to establish the 
equivalence of H and H. Consider the one- dimensional Schrodinger eigenvalue problem 

h 2 

ip"(x) - gx 4 ip(x) = Eip(x) (4) 

2m 

associated with the non-Hermitian Hamiltonian (|T)). where the boundary conditions on ip{x) 
are that as lim| 9; |_ >oo '0(a:) = if — | < args < and — n < argx < — These boundary 
conditions do not include the real-x axis and they require that the differential equation (@J) 
be solved along a contour whose ends lie in the above wedges in the complex- a; plane (see 
Fig.H}. We will use here the same complex contour that Jones and Mateo employed in their 
operator analysis of the Hamiltonian (0) pi 21]: 



x 



-2iLy/l + iy/L, (5) 



where y runs from — oo to oo along the real axis. This contour is acceptable because as 
y — > ±oo, argx approaches —45° and —135°, so the contour lies inside the Stokes' wedges. 
In (JSJ) L is an arbitrary positive constant having dimensions of length. In terms of the 
parameters of H in (0) the fundamental unit of length is [h 2 / (mg)] 1 ^ 6 . Thus, 

/ h 2 \ 1/6 

L = X ( — — j , (6) 
\mgj 

where A is an arbitrary positive dimensionless constant. When we change the independent 
variable in (jij) from x to y according to (0), the Schrodinger equation (@J) becomes 

~H 1 + t) r(v) ~ tL* (v) ~ lt,3L ' ( 1 + t) 2 ^ = my) - (7) 

Next, we perform a Fourier transform of ((Zj). We define 

poo 

f(p)= / dye-^ n f(y), (8) 
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so that the Fourier transform of f'(y) is ipf(p)/h and the Fourier transform of yf(y) is 
ihf'(p). Then, the Fourier transform of the Schrodinger equation (J2J) is 

Expanding and simplifying this equation, we obtain 
- 1Q 9 L^4>"( P ) + f-M + 32gL 3 h)j 4>'(p) +(£-^L- 16 9 lA fo) = E$(p)(10) 

[Note that the variable p used here is not the same as the variable p used in (pQ). Here, 
considered as an operator, p represents —ih-^, whereas in ((TJ) p represents —ih-^.} 

This equation is not a Schrodinger equation because there is a one-derivative term. How- 
ever, we can eliminate this term by performing a simple transformation: 

0(p) = e Q(p)/2 $(p). (11) 

The condition on Q(p) for which the equation satisfied by $(p) has no one-derivative term 
is a first-order differential equation whose solution is 

2L 1 

QM -T p -96^m p3 - (12) 



It is interesting that is precisely the operator found in Ref. [iSf 1 Substituting this 
expression for Q gives the Schrodinger equation satisfied by $(p): 

Finally, we make the scaling substitution 



p = zL v /32mg, (14) 

to replace the p variable, which has units of momentum, by z, which is a coordinate variable 
having units of length. The resulting eigenvalue equation, posed on the real- ,2 axis, is 

- — $"(*) + \-h\j^-z + Agz A \ $(z) =E${z). (15) 
2m \ V m I 

We emphasize that while z has dimensions of length, it is not a conventional coordinate 
variable because it is odd under the discrete transformation of time reversal. 

Observe that the eigenvalue problem (|15|) is similar in structure to that in (J3J). [Equation 
(115(1 is not dual to (jlj because it is still weakly coupled.] However, the potential has acquired 



Note that e^ p ' is not the CP operator that could, in principle, be used in a similarity transformation to 
produce the Hermitian Hamiltonian from the non-Hermitian 'PT-symmetric Hamiltonian Ijljl. One can 
only use to transform the non-Hermitian Hamiltonian H to Hermitian form if H has first been 

written in terms of the real variable y. 



5 



a linear term, and since this linear term is proportional to h, we may regard this term as a 
quantum anomaly. The linear term has no classical analog because the classical equations 
of motion are parity symmetric. The breaking of parity symmetry occurs at large values of 
x where the boundary conditions on the wave function ip(x) are imposed. Because we have 
taken a Fourier transform to obtain the Schrodinger equation (|15|). this parity anomaly now 
manifests itself at small values of z. 

The Hamiltonian H for which (|15|) is the eigenvalue problem is 

H=^-h\[^z + 4gz\ (16) 
2m V tti 

This Hamiltonian is Hermitian in the Dirac sense and is bounded below on the real-z axis. 
Furthermore, it is also PT-symmetric. This is because at every stage in the sequence of 
transformations above, VT symmetry is preserved. However, while z and p are canonically 
conjugate operators satisfying [z,p\ = i, the new variable z behaves like a momentum rather 
than a coordinate variable because z changes sign under time reversal. 



III. BOUND STATES A DIRECT PHYSICAL CONSEQUENCE OF THE 
ANOMALY 

If we generalize the Hamiltonian to include a harmonic term in the potential, 

H = —p 2 + ^-x 2 - gx\ (17) 

then the same differential-equation analysis used in Sec. [TTI straightforwardly yields the fol- 
lowing equivalent Hermitian Hamiltonian: 

S = f -J*, + 4,(z--f)'. (is) 

2m \ m \ 8g J 

This result is given in Refs. 0, . Observe that for these more general Hamiltonians the 
form of the linear ano maly term remains unchanged from that in (fT3j) . 

In an earlier paper 14] it was shown that the Hamiltonian (fT7j) exhibits bounds states. 
In a particle physics model a bound state is defined as a state having a negative binding 
energy. In the context of a quantum-mechanical model we define bound states as follows: 
Let the energy levels of the Hamiltonian be E n (n = 0, 1, 2, . . .). The renormalized mass is 
the mass gap; that is, M = E\ — E$. The higher excitations must also be measured relative 
to the vacuum energy: E n — E {n = 2, 3, 4, . . .). We say that the nth higher excitation is a 
bound state if the binding energy 

B n = E n -E -nM (19) 

is negative. If B n is positive, then we regard the state as unbound because this state can 
decay into n 1-particle states of mass M in the presence of an external field. 

It was observed numerically in Ref. 0] that for small positive values of g the first few 
states of H in ffTTj) are bound. As g increases, the number of bound states decreases until, 
when g//i 3 is larger than the critical value 0.0465, there are no bound states. In this paper 
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a rough heuristic argument was given to explain why there is such a critical value. This 
argument is difficult to formulate because the non-Hermitian Hamiltonian is evaluated for 
x in the complex plane. When x is complex, one cannot use order relationships such as > 
or <, which only apply to real numbers. 

However, now that we have established that H in (|17|) has the same spectrum as the 
Hermitian Hamiltonian in (|18|) . it is easy to understand the appearance of bound states. 
Furthermore, as we will now show, the bound states are a direct consequence of the linear 
anomaly term. To probe the influence of the anomaly, let us generalize (fTHj) by inserting a 
parameter e that measures the strength of the anomaly term: 



H 



P 2 



ehy2g z + Ag 



1 

&9 



(20) 



where for simplicity we have set m — \i = % = 1. 

If we take e = 0, then there is no anomaly term and the potential is a symmetric double 
well. The mass gap for a double well is exponentially small because it is a result of the 
tunneling between the wells. Thus, the renormalized mass M is very small. Therefore, 
B n in (jTUj) is positive and there are no bound states. In Fig. |2] we display the double- well 
potential and the first several states of the system for the case g = 0.046 and e = 0. Note 
that the lowest two states have a very small splitting. 

If e = 1, the double- well potential is asymmetric and the lowest two states are not 
approximately degenerate. As a result, bound states can occur near the bottom of the 
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FIG. 2: Potential of the Hermitian Hamiltonian (|2Uj) plotted as a function of the real variable z for 
the case e = and g = 0.046. The energy levels are indicated by horizontal lines. Because e = 
there is no anomaly, and the double-well potential is symmetric. Therefore, the mass gap is very 
small and, as a consequence, there are no bound states at all. The occurrence of bound states can 
be attributed to the anomaly. 
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potential well. The higher-energy states eventually become unbound because, as we know 
from WKB theory, in a quartic well the nth energy level grows like n 4//3 for large n. As g 
becomes large, the number of bound states becomes smaller because the depth of the double 
well decreases. For sufficiently large g there are no bound states. In Fig. |H] we display the 
potential for e = 1 for g = 0.046. For this value of g there is only one bound state. 

To display the bound states we simply plot the value of the binding energy B n as a function 
of n. For example, in Fig. HJwe display the bound states for e = 1 and g = 0.008333. Note 
that for these values there are 23 bound states. Observe also that the binding energy B n is 
a smooth function of n. 

It is interesting that the bound-state spectrum depends sensitively on the anomaly term 
in the Hamiltonian (j2*Uj) . If e is slightly less than 1, the first few states become unbound, as 
is shown in Fig. In this figure g = 0.008333 and e = 0.9. 

If e is slightly greater than 1, the binding energy B n is not a smooth function of n for 
small n. In Fig. El we plot B n as a function of n for g = 0.008333 and e = 1.1. Note 
that for these values of the parameters there are 30 bound states. Figures 01 El and El are 
strikingly different, and demonstrate the extreme sensitivity of the bound-state spectrum to 
the anomaly term. 

IV. PATH-INTEGRAL DERIVATION OF THE PARITY ANOMALY 

A path integral is used to represent the partition function of a quantum field theory. A 
path integral is much more complicated than a Schrodinger equation, and thus the elemen- 
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FIG. 3: Asymmetric potential well plotted as a function of the real variable z for the Hermitian 
Hamiltonian (|20[) with e = 1 and g = 0.046. The energy levels are indicated by horizontal lines. 
There is one bound state. 
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tary differential-equation methods applied in Sec. HP cannot be used directly for a quantum 
field theory To prepare for our study of a PT-symmetric — 4 field theory in Sec.|3 we need 
to establish the equivalence of H in (Q) and H in ([To*)) by using path-integration methods. 
For the Hamiltonian the Euclidean functional integral for the partition function Z is 



Dcf) exp 



c 



dt 



m 



(21) 



where the normalization factor is understood. A functional integral is an infinite product of 
ordinary integrals, one for each lattice point in the discretized representation. The g4> A term 
in the exponent would cause each of these integrals to diverge if the integration path lay on 
the real axis. To make these integrals converge, the contour of integration at each lattice 
point must approach infinity inside a pair of wedges having an angular opening of 45° and 
centered about the angles —45° and —135° in the lower-half complex plane. The subscript 
C on the functional integral sign indicates that the path integral is taken along a complex 
contour. 

Our goal is to transform the functional integral (}2"Tj) into a conventional functional integral 
in which the contour of integration runs along the real axis rather than in the complex plane. 
We expect to find that the action in the exponent will correspond to H in 

Our approach is as follows: First, in Subsec. IIV Al we transform the continuum func- 
tional integral in (}2"T]) directly using (jSj). We discover that if we proceed formally without 
recognizing that a functional integral is a subtle construct involving the limit of a discrete 
product of ordinary integrals, the resulting equivalent Hermitian Hamiltonian that we ob- 
tain does not contain the linear anomaly term in H in (fTrjj) . Thus, this derivation produces 
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FIG. 4: Binding energies B n = E n — Eq — nM plotted as a function of n for g = 0.008333 and 
e = 1. A negative value of B n indicates a bound state. Observe that there are 23 bound states for 
these parameter values. Note that B n is a smooth function of n. 



9 



the classical (h — > 0) limit of H. Hence, this derivation is equivalent to the geometrical- 
optics approximation (zeroth-order WKB). To obtain the anomaly, we must discretize the 
functional integral, and this discretization requires great care. We explain and motivate the 
discretization procedure in Subsec. II V Bl and carry it out in Subsec. II V CI 



A. Formal Derivation Correct to Order h 1 



o 



We begin by making a substitution analogous to that in (JHJ), 



(f>(t) = -2iL^/l + i^(t)/L. 
This substitution introduces a functional Jacobian in the form of a square root: 

Dtp 

D0 = . = 

Det y/l + ii{>/L 

The new functional integral over the ip variable is 



GXO \ / dit 

Det y/1 + iip/L \ h 



m ^ 2 (t) 
~2 1 + iip(t)/L 



- 16gL A [l + i 



L 



(22) 



(23) 



• (24) 
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FIG. 5: Binding energies B n plotted as a function of n for g = 0.008333 and e = 0.9. Note that the 
first five states have now become unbound and B n is not a smooth function of n for n < 6. The 
next twelve states are bound, and in this region B n is a smooth function of n. A comparison of 
this figure with Fig. 0] shows that the bound-state spectrum is exquisitely sensitive to the strength 
of the linear anomaly term. 
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This is a conventional functional integral in the sense that the field ip(t) is real and the path 
of integration lies on the real axis rather than in the complex plane. 

Next, we exploit the functional- integral version of the simple integral identity 



i = -L r dte~\ A2 ^ 2 (ReA 2 >0) 



A 

and obtain (apart from an overall normalization factor) 



(25) 



Dip I Dn exp {-- / dt 



l + iil>(t)/L ( 



2m 



\*{t) - B(t) 



m ip 2 (t) 
'~2 1 + iip(t)/L 



-!QgL 4 (l + iij(t)/L 



(26) 



This integral converges because Re[l + itp(t)/L] = 1, which is positive. Both functional 
integrals in ()26|) are conventional real integrals. 

We eliminate the term containing the square of ip(t) by choosing 



B(t) 



With this choice the exponential in 



imip(t) 
~~ l + ty(t)/L' 

simplifies to: 



(27) 



\ 1 



— 7r 2 (t)(l + #(t)/L) - iir(t)i>(t) - lQgL A (l + iiP{t)/L Z 



(28) 
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FIG. 6: Binding energies B n plotted as a function of n for g = 0.008333 and e = 1.1. Note that 
there are 30 bound states and that B n is not a smooth function of n when n is small. 
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Next, we integrate by parts in order to transfer the derivative from ip(t) to 7r(t) and in 
doing so discard the surface term: 



dt7p(t)n(t) = / dtn(t)ij}(t). 



(29) 



After interchanging orders of integration and completing the square, we get 

w .7 i 



Dti exp 
x / Dip exp 



11 J dt \Q4gL 2 



tV 2 



(t) + 



' vr 4 (t) - Lk{t) + ' 



256#m 2 L 4 



QAgmL 



16gL 2 



dt %fi{t) -iL + 



64gmL 3 



At) + 



32gL 2 



*(*) 



(30) 



We integrate the n(t) and 7r 2 (t)7r(t) terms and ignore the surface contributions. 
Finally, we rescale 7r and ifi in analogy with (|14j) : 



7T 



(t) = Ly/32mg<p(t), 



L^32mg 

and evaluate the Gaussian integral over p. The result is 



*>(*), 



.D^ 9 ex P 



-iy 'dt(^(i)+w(()) 



(31) 



(32) 



which is the Euclidean functional integral for the classical (anomaly-free) version of the 
Hamiltonian H in flloj) . 



B. Discretization of the Function Integral 

Evidently, the analysis in Subsec. II V Al is not delicate enough to recover the linear anomaly 
term of H in (|lb)). In order to obtain this anomaly term it is necessary to discretize the 
functional integral (J21j) . In order to discretize, we replace <f>(t) by the usual lattice expression 
(4> n +i — (fin)/®, where a is the lattice spacing. Furthermore, we must replace <j>(t) by the 
average \{(fi n+ \ + (fin) to preserve the time independence of the equal-time commutation 
relation in the underlying quantum-mechanical theory. Our purpose in this subsection is to 
explain the reason behind this latter substitution. 

Consider a one- dimensional quantum-mechanical Hamiltonian of the form 

H=\p 2 + V(q). (33) 

The Heisenberg equations of motion of the operators p(t) and q(t) have the same form as 
the classical equations in 0: 

x{t)=p, p{t) = -V'[q{t)}. (34) 
To discretize the differential equations in f)34j) . we write 

%n+l -En Pn+1 Pn Pn+1 Pn t rf I "^n+1 \ /nr\ 

—a = ^^' -^— = - V { 2 J' (35) 
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where a is the lattice spacing. 

One can verify that (|35|) is correct by following the procedure introduced in Ref. (lif . 
The right side of the first equation in (|3*5|) is a function of p n +\ + Pn, so the commutator of 
this equation with p n +i + p n is 

[X n+1 - X n , Pn+l + Pn] = 0. (36) 

Also, the right side of second equation in (|33j) is a function of x n +i + x n . Thus, commuting 
x n+ \ + x n with this equation gives 

[X n+ i + X n , Pn+l ~ Pn] = 0. (37) 

Adding (|3T)j) and (|3*7|) then gives 

[x n+1 ,p„+i] = [x n ,p n ] . (38) 

This establishes the crucial result that the equal-time commutator is exactly preserved in 
time. Thus, the discretization scheme in (|3*5)l is exactly unitary. If a discretization scheme 
other than that in (|35|) had been used, there would be a small violation of unitarity. Ev- 
idently, to avoid violating unitarity on the lattice it is essential to replace a local function 
of x(t), say f[x(t)} by /[| (x n+1 + x n )] rather than by f(x n ). Of course, the violation of 
unitarity vanishes in the continuum limit a — > 0, but we need to preserve unitarity in the 
lattice version of the theory. 

To discretize an action, we follow the same approach as that used to obtain (|35|) . Thus, 
for the Lagrangian L = \x 2 (t) — V[x(t)\ we discretize the action J dt L as follows: 



/VI v J 



(39) 



To verify that this is the correct way to discretize the action, we vary this discrete action 
with respect to x n to obtain the lattice equations of motion: 

1 Qj Ct 

- (2x n - X n+ i - X n -i) - -V [\ (x n+ i + x n )] - -V [\ (x n + X n - X )\ = 0. (40) 
To solve this equation, we substitute 

Xn+l ~ x n Pn+l + Pn 



(41) 



a 2 

which is the first equation in (J35|) . This reduces (|4U|) to the form 

b n + b n _x = 0, (42) 

where 

K = P n + l-Pn +v ,f x n+ x+X n \ (43) 



The solution to (J42D is 

b n = c(-l) n , (44) 
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FIG. 7: Stokes' wedges in the lower-half complex-c/> n plane for each of the one-dimensional complex 
integrals used in (|45[). The integrands of these complex integrals vanish exponentially fast as 
\(j) n \ — ► oo within these wedges. Also shown is the complex contour in Q46J1 . 

where c is an arbitrary constant. However, if c is nonzero, then in the continuum limit where 
the lattice spacing tends to the solution in ()44|) becomes infinitely oscillatory. To avoid 
the appearance of an infinite-energy solution, we must require that c = 0. We therefore 
obtain the result that b n = 0, and we have recovered the second equation in We have 
reproduced both equations in (|Ho]l and hence conclude that the discretization scheme used 
in is correct. The lattice average used on the right sides of (f3T)|) and in the second term 
of ()39|) is completely consistent with Moyal ordering [l6| . 



C. Derivation of the Parity Anomaly 

Deriving the Hamiltonian H in from the partition function Z in (|21j) requires great 
care. First, we must discretize the functional integral (J2Tj) representing Z by following the 
procedure in Subsec. IIVBI 

Z = Yl / #n exp |-| ^ (4> n+1 - (j) n ) 2 - ^ (<f> n+1 + n ) 4 I , (45) 



c 



where a is the lattice spacing. For each of the one-dimensional integrals in ()45j). the contour 
C begins in a wedge in the lower complex plane centered about —135° and terminates in 
the PT-symmetric wedge centered about —45° (see Fig. EJ). 

For each of the complex integrals in (|45|1 we make a substitution analogous to that in (j3J): 

<P n = -2iLy/l + iif> n /L, (46) 

where L is an arbitrary positive constant having dimensions of length. At each lattice point 
this substitution introduces a Jacobian in the functional integral; 

ii dlp n 

d<Pn = /a -L ' / IT ' (4?) 
A/1 + Ityn/L 
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Thus, an exact transcription of (|4*5j) is 



n 



where we use the shorthand notation 



exp 



a ( m ipl 



16gL*A n 



2AI 



(48) 



1pn = ~ (VVH-1 _ 

a 



and 



A n = - (yJ\+ilj) n+1 /L + a/1 + ilpn/L^J . 



(49) 



(50) 



Note that in (j48j) the path of integration of each of the integrals lies on the real axis. 

We now use the integral identity (J2HJ) to introduce the conjugate auxiliary field 7i n at each 
lattice site n. To apply this identity we choose 



A 1 



— A 2 
mn 



t = 7T„ 



B = B r , 



(51) 



where B n is as yet unspecified. The result is 



dip„ 



dlTr 



A. 



x exp 



a 

n 



Mm a/I + itj) n /L 



(52) 



Making no approximations, we simplify the Jacobian in (J52|) as follows: 

A n 



a/1 + il/jn+l/L + a/1 + #n/£ 



a/1 + # n /£ 



2 a/1 +ii/> n /L 



1 a/1 + iljjn+l/L - a/1 + 

2 a/1 + #„/L 

z(^ n+ l - V»n) 



+ 
1 + 



1 + a 



4LA nV /l + # n /L 



Next, we choose B n : 



B„ 



imipr, 



(53) 



(54) 



[Note that the choice of sign on the right side of (JSlj) is arbitrary. If we replace i by —i in 
this equation, then at the end of the calculation the sign of the linear term in the Hermitian 
Hamiltonian will be reversed. Reversing this sign has no effect on the energy levels of the 
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Hamiltonian.] With this choice, we can approximate for small lattice spacing a the Jacobian 
factor in the last line of (|53|): 



1 + a- 



iip ri 



ALA n ^l+i^ n /L 
Next, we use the identity 

aAl 



l + ^ + 0(a 2 ). 
AmL v ' 



dn n 7i n exp 



2mh 



{7l n - B n ) 



dix n B n exp 



nA 2 

[7T„ - ±f„) 



2mh 



(55) 



(56) 



which holds because the exponent is an even function of 7r n — B n . This identity implies that 
(|52|) with the Jacobian in (jo3j) can be rewritten as 



n 



dip. 



2nTim 



l + ^+0(a 2 ) 
AmL v ' 



x exp 



i^n^n - lQgL 4 A 4 



h \ 2m 



(57) 



Next we promote the Jacobian factor to the exponent by noting that for any number w, 
we have 1 + aw = e aw with an error of order a 2 . We also simplify terms in the exponent by 
ignoring terms of order a 2 . The result is 



n 



dlpr, 



dir. 



" w 2nhm ^ I h 



1 + ilp n /L) — ^ + 27T n ^ r , 

2m 



+ O a) 



■l6gL*(l + ii; n /Ly 



AmL 



(58) 



where we have performed a summation by parts to obtain the i7r n ip n term. 
Next, we perform a scaling like that in (|31|): 



7r n = ^/32mgLip n , 



-Pn, 



(59) 



^/32rngL 

and continue as we did in Subsec. II V Al After completing the square and integrating over 
p n , we obtain the final result for the lattice version of the continuum functional integral: 



Dip exp 



(60) 



This is precisely the functional integral for the theory described by the Hermitian Hamilto- 
nian H in (|16|). 

V. GENERALIZATION TO QUANTUM FIELD THEORY 

The Euclidean functional integral for the partition function Z of a d-dimensional — <f) 4 
quantum field theory is 



Dtp exp 



c 



d d x{\{V<P) 2 -g^) 



(61) 
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where the normalization factor is understood and we have adopted natural units where 
h = 1. As in the one-dimensional quantum field theory in (|21|). the contour of integration 
at each space-time point must approach infinity inside a pair of wedges having an angular 
opening of 45° and centered about the angles —45° and —135° in the lower-half complex 
plane. The subscript C on the functional integral sign indicates that the path integral is 
taken along a complex contour. 

As in the quantum-mechanical case, our goal is to transform the functional integral (|61j) 
into a conventional functional integral in which the contour of integration runs along the 
real axis rather than in the complex plane. Following the procedures in Sec. IIVI we make 
the substitution 



<j>{x) = -2iy/l+iip(x), (62) 

where we have set L — 1. This substitution introduces a functional Jacobian in the form of 
a square root: 

Dcp = ^^t-^j . (63) 
Det v 1 + tip 



The new functional integral over the ip variable is 



Dip 



exp s — d d x 



Det y/1 + iip [J 1.2(1 + z^>) 



(WO 



2 



- 16#(1 + iip) 



(64) 



This is a conventional functional integral in the sense that the field ip(x) is real and the path 
of integration lies on the real axis rather than in the complex plane. 

From the form of the path integral in (|64j) we can see that even though the complex 
number i appears in the integrand, Z is real because the change of variable ip — > — ip 
changes the sign of i. Thus, the ground-state energy density is real. 

This same argument implies that the 2n-point Green's functions are real and the (2n+ 1)- 
point Green's functions are imaginary. The n-point Green's function G n (xi,X2,xs, . . . ,x n ) 
is constructed from the moments of products of (p fields: 



(O|0(xi)0(x 2 )0(x 3 ) • • -(p(x n )\0) = I D(p(p(x 1 )(p(x2)(p(x 3 )---(p(x n ) 

d d x{\{Vcp) 2 -g^) 



c 

x exp 



(65) 



After making the substitution ()62j). we can apply the above symmetry argument to show 
that G n is real when n is even and imaginary when n is odd. 

Even though we have shown that the partition function Z in ()64j) is real, we are still 
unable to identify a Hermitian field-theoretic Lagrangian that is equivalent to the — <p A non- 
Hermitian Lagrangian. This is because the functional integral ()64|) has a Jacobian factor, 
and this factor is complex. By following the procedure in Sec. IIVI we were able to eliminate 
this Jacobian factor for the case of a one-dimensional quantum field theory. It is not obvious 
how to proceed in the case of higher dimensions. We now propose three possible approaches 
to solving this problem. 

1. The simplest approach is to follow Sec. IIVI and introduce a new field 7r(x) by using the 
identity (|23j) . The advantage of this procedure is that, at least at the classical level 
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where we do not perform the point-splitting described in Subsec. IIV CI the Jacobian 
factor in the functional integral cancels. We obtain the (i-dimensional analog of ()26|) : 



Dip Dn exp<- d d x 



There are now several possibilities for the choice of B. The choice 



.(66) 



[) l + iiP{x) [ ' 

eliminates the [Vi?/>(£)] 2 term from (}66|) but the resulting action is nonpolynomial in 
ip and thus the integral over ip cannot be evaluated analytically. Furthermore, the 
resulting action is not manifestly Lorentz covariant. The more symmetric choice 

= (68) 
v ; l + iip(x) K ' 

eliminates the (V^) 2 /(l + iip) term completely, but again the resulting action is non- 
polynomial and it is impossible to perform the ip integral. 

2. A second approach is to introduce the (i-dimensional vector field 7r M (// = 1, . . . d) . Now, 
using the identity (J25j) d times, we obtain an action that is a quadratic polynomial in 
ip. However, we still cannot perform the ip integration because the Jacobian factor no 
longer cancels. Indeed, there are now d — 1 functional determinant factors of a/1 + iip 
in the numerator. Furthermore, had we kept the scale factor L in the change of vari- 
ables the resulting Jacobian would have explicit L dependence, even though the 
partition function is clearly independent of L. Therefore, this integral representation 
of the partition function is clumsy and is likely to be intractable. 

3. We believe that the most promising approach is to introduce a single scalar field p(x) 
by formally rewriting Z in (J64j) as 



Det v/V M (l + #)V M 

Dip . 

Det y/1 + iip 

x J Dp exp j- J d d x [|(Vp) 2 (l+#)-iV M pV^-16^(l + #) 2 ]|-(69) 

The advantage of this representation of Z is that the action is quadratic in ip and is 
Lorentz covariant. The Jacobian remains complicated, but it may become analytically 
tractable if we introduce a Grassmann integration variable 77 to promote the functional 
determinant to the action. After integrating out the ip field, the action will now depend 
on the scalar boson field p and the fermion field rj. The resulting Hermitian action will 
bear a strong resemblance to supersymmetric actions. This is not surprising because, 
even in the one-dimensional case, the equivalent Hermitian Hamiltonian H in |T6|) 
strongly resembles a supersymmetric quantum theory. In a future paper we hope to 
present a thorough analysis of this quantum theory. 
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